Squeezed states and quantum chaos 
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We examine the dynamics of a wave packet that initially corresponds to a coherent state in the 
model of quantum kicked rotator. This main model of quantum chaos, which allows for a transition 
from regular to to chaotic behavior in the classical limit, may be realized in experiments with cold 
atoms. We study the generation of squeezed states in the quasiclassical limit and in a time interval 
when quantum-classical correspondence is yet well-defined. We find that the degree of squeezing 
depends on the degree of local instability in the system and increases with the Chirikov parameter 
of stochasticity. We discuss the dependence of the degree of squeezing on the initial width of the 
packet, the problems of stability and observability of squeezed states at the transition to quantum 
chaos, as well as the dynamics of wave packet destruction. 



I. INTRODUCTION 

At present the problem of squeezed quantum states 
generation draws a lot of attention, both from the 
standpoints of pure knowledge and possible applications 
10, H, ||]. Most often the topic is squeezed states of the 
electromagnetic field. If in the simplest case we consider 
a single-mode quantum field, which is described by the 
creation and annihilation a operators, the variances of 
quadrature field operators oi = a-\-a^ and 02 = ~i{a—a^) 
satisfy the uncertainty relation AaiAa2 > 1, where the 
equality holds for a coherent state or vacuum. Then, in 
these simple terms, a squeezed state is a state for which 
the variance of one of the quadrature components is less 
than unity. Quantum fluctuations, determined by the un- 
certainty relation, can be represented diagrammatically 
in the oi — 02 plane by a circle for a coherent state or 
by an ellipse for a squeezed state. In a more system- 
atic description of squeezing, the quantum-noise ellipse 
is determined in terms of the projection onto the ai — 02 
plane of the horizontal section of the Wigner distribu- 
tion function, which gives the quasiprobability distribu- 
tion for measuring the quadrature field components ^ . 

A typical situation in experiments on generation of 
squeezed states is one in which a large number of pho- 
tons participate in a nonlinear interaction and the am- 
plitude of quantum fluctuations is small compared to the 
averaged values of the observables In this case 

the common approach in explanation of squeezing is to 
use the semiclassical setting, where the Wigner quantum 
function is actually replaced by a classical distribution 
function and instead of examining the dynamics of the 
quantum-noise ellipse one considers the evolution of the 
classical phase volume ^, ^ . 

For quite a long time it has been known that squeez- 
ing is enhanced if a system is close to a bifurcation point 
between two different dynamical regimes ^ ||, ^ ■ The 
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increase of squeezing in such conditions was considered, 
e.g., for the parametric interaction of light waves § and 
for the interaction of Rydberg atoms with an electromag- 
netic mode in a high-Q cavity in a dynamical regime close 
to the separatrix |^. The following simple argument 
is used to explain the increase of squeezing near a bi- 
furcation point. Quantum fluctuations build up for the 
physical variable that is unstable near the threshold. As 
a result there is nothing to stop the strong squeezing for 
the conjugate variable because of conservation of a phase 
volume in a nondissipative system 

It must be noted at this point that a number of re- 
searchers (see Refs. (|, @, ||) studied the increase 
of squeezing near the instability threshold in different 
optical systems with regular dynamics. However, it is 
well known that strong (exponential) deformation of the 
phase volume is one of the main manifestations of dynam- 
ical chaos in classical systems [Q| . The physical reason for 
such strong deformations of the phase volume is the local 
instability of chaotic motion, which usually manifests it- 
self within a wide range of values of the control parameter 
of dynamical system. According to the correspondence 
principle, in the quasiclassical limit a quantum system 
must mimic the properties of a classical system. Thus, 
it is quite natural to expect the increase of squeezing at 
the transition to quantum chaos. On the other hand, 
in a quantum mechanical description we speak only of 
the dynamics of wave packets, whose center moves al- 
most along a classical trajectory in the course of a cer- 
tain time interval. Hence, the strong deformations of the 
phase volume, which accompany the transition to chaos, 
must manifest themselves in squeezing along a certain 
directions in phase space. 

As far as we know, the generation of squeezed states 
in systems with chaotic dynamics was first examined 
in Refs. J&_P|j 0- By employing the 1/A^-expansion 
method p,^! Ij (here N is the number of quantum states 
participating in the dynamics of the system) it was found 
in Refs. || that the squeezing of light increases signif- 
icantly at the transition to chaos during the time inter- 
val for which quantum-classical correspondence is well- 
defined Q. This result was first obtained in Refs. |8[ ||] 
for the generalized A'^-atoms Janes-Cummings model, 
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which allows a transition from regular to chaotic dynam- 
ics in the limit N ^ oo . Later the increase of squeez- 
ing at transition to chaos was found for arbitrary single 
mode quantum optical systems p^ . The squeezing of 
wave packets in conditions of quantum chaos was also 
briefly discussed in However, the main results of 

works §, g were obtained by using a form of pertur- 
bation theory (the 1/A^-expansion) . Therefore it is im- 
portant to investigate the generation of squeezed states 
at transition to quantum chaos in a simple quantum sys- 
tem whose time evolution can be found explicitly. 

In the present paper we study the generation of 
squeezed states in the time evolution of an initially Gaus- 
sian wave packet in the model of quantum kicked rotator. 
This model was first introduced by Casati et al. and 
at present is the main model in studies of quantum chaos 
(for review see |l^, Q ) . The quantum rotator model 
is attractive mainly for two reasons: (i) the classical limit 
for this model is a well-known standard map jiot , and (ii) 
it is fairly easy to study numerically the dynamics of the 
model for a large number of quantum levels. 

We examine the dynamics of narrow Gaussian pack- 
ets in a rotator with 2^^(~ 10^) levels. We define 
squeezing for the generalized quadrature operator Xg — 
aexp{~i9) + e'xjp{9), where is a real parameter. It 
is just Xg that is observed in the homodyne detecting 
scheme, where 6 is determined by the phase of the ref- 
erence beam Q. We will see that as long as the wave 
packet is localized, the degree of squeezing correlates 
well with the degree of local instability in the system. 
Here the greater the instability, the stronger the squeez- 
ing achieved in a shorter time interval. Squeezing is much 
stronger for quantum chaos than it is for regular motion. 
We will also see that the narrower the initial wave packet, 
the higher the degree of squeezing that can be achieved. 
We attribute this to the fact that a narrow wave packet 
is closer in its time evolution to the classical trajectory 
than a broad one, with the result that it is more sensi- 
tive to local instabilities in the motion, which leads to 
stronger squeezing. 

We will also consider the problem of stability and ob- 
servability of squeezed states in the transition to chaos. 
More precisely, we will study the time dependence of the 
optimum values of the phase of the generalized quadra- 
ture operator Xg for which the squeezing is maximal (so- 
called principal squeezing ||2^, ^). We will show that 
in conditions of strong chaos and for long enough time, 
the optimum values of the phases change dramatically 
even under a small perturbation of the parameters of the 
initial Gaussian wave packet. Such a squeezing regime 
is unstable and difficult to observe. On the other hand, 
our results suggest that for weak chaos squeezing is fairly 
stable. 

We will also briefly discuss the dynamics of disintegra- 
tion of wave packets at chaos. Here we will show that a 
typical scenario of disintegration of an initially localized 
wave packet in chaos consists of two stages: the initial 
spread of the wave packet, and the catastrophic disinte- 




gration of the packet into many small subpackets. Here 
our results agree on the whole with the results of Casati 
and Chirikov [g. 

Note that earlier the dynamics of narrow Gaussian 
wave packets in the quasiclassical limit was studied nu- 
merically for the models of a quantum kicked rotator 
quantum kicked top and quantum Arnold 
in connection with the problem of quantum- 
classical correspondence at quantum chaos. However, in 
these papers the generation of squeezed states was not 
considered. 

The model of a quantum kicked rotator is very popular 
in theoretical studies of quantum chaos. On the other 
hand, recently possibilities of implementing variants of 
this model in optical systems have been discussed [ p6[ . 
Moreover, the model of quantum rotator has been real- 
ized in experiments on the interaction of laser beam and 
cooled atoms Hence our results on the increase of 

squeezing at the transition to quantum chaos in a kicked 
quantum rotator are also directly related to experimen- 
tally realizable systems. 

The plan of the paper is as follows. In Sec. II we dis- 
cuss the quantum standard map and find how to calculate 
principal squeezing. The method used in numerical sim- 
ulations is developed in Sec. [H, and the main results on 

Finally, 



the dynamics of squeezing are given in Sec. IV 
m Sec. we draw the main conclusions and consider the 
possibility of verifying our results in experiments. 



II. THE MODEL OF QUANTUM KICKED 
ROTATOR AND SQUEEZED STATES 

Consider the model of a quantum rotator with periodic 
delta-kicks. Here we follow the notation of Ref. |2^. The 
Hamiltonian is 



H 



P 



5p{t/T)mL ujqCOs{x), 

+ 00 



2mL2 



(1) 



where x is the cyclic variable with a period 27r, L is the 
characteristic size of the rotator, m is the rotator mass, 
and uiQ is the frequency of linear oscillations. The func- 
tion Sp{t/T) describes a periodic sequence of kicks with 
a period T, where 6{x) is the Dirac (5-function. Let us 
introduce new variables 



T 



(2) 



and measure time in units of T, i.e., t — s- t/T. Then the 
Schrodinger equation takes the form 



ih 



(9* 
'dt 



6p{t)a cos (x)\l/. 



(3) 
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Due to the periodicity of ^{x) in x the solution of Eq. 
can be written as foUows 



Akit) = 



1 



+ 00 



1 



/c— — 00 
2ir 



^{x)e~"'''dx. 



(4) 



Using the standard procedure 
quantum map in the form 



we obtain the 



Up — exp 



(5) 



i(3 , 
2^^ 



Ux — exp ( — cos [x) 



where is the value of the wave function at the time 
moment immediately after the rt-th kick. The time evo- 
lution of the wave function in the map (|^) is determined 
solely by two parameters, a/Ti and j3fi. Since Up is di- 
agonalized in the p-representation, Ux is diagonalized in 
the a;-representation, and the transition between x— and 
^-representations is given by the Fourier transformation 
(0), the map (0) actually reduces to 



n+l 



(x) = UxF-^UpF-^r.{x), 



(6) 



where F and F^^ are the direct and inverse Fourier trans- 
forms. 

Sometimes it proves useful to use the quantum map 
written in terms of the probability amplitudes A^ of tran- 
sitions between the unperturbed levels of the rotator . 
Combining Eqs. (||) and we obtain 



A 



(n+l) 



E 



Fk Ai 



(7) 



Fk 



k — m 



exp 



where Ji{z) is the Bessel function of order I and argument 
z, and the superscript (n) on the variable A stands for 
the number of the kick. Bearing in mind that the Bessel 
functions with | ? |> 2: rapidly decrease with increasing 
I, we see from Eq. (|^) that, with exponential accuracy, 
in the course of a single kick 2aTi~^ unperturbed rotator 
levels are captured. Below we consider the case where 
a/?i is large, which is typical of quantum chaos problems. 

In the classical limit the Hamiltonian (|^) reduces to 
the standard map 



Pn+l = Pn ~ KsmXn+l 



Xn+l 



n + Pn (mod 27r), 
(8) 

where P„ = Ppn, the subscript n denoting the values of 
X and P immediately after the n-th kick, and K = afi is 
the Chirikov parameter^. Strong and global chaos sets in 



The unusual form of the standard map (p|) is due to the fact that 




FIG. 1: Phase portrait of the classical standard map for K ■ 
0.8. 



for K > 1. For K < 1 the larger part of the phase plane 
is filled with regular trajectories, although small regions 
with local chaos exist no matter how small K may be . 
The phase portrait for the map (^) at iiT = 0.8 is depicted 
in Fig. 0. The chaotic layer lies near the separatrix of 
the main resonance, which passes through the hyperbolic 
points (±7r, 0). In our calculations we usually take a wave 
packet whose center is located near a hyperbolic point 
but yet in a region of regular motion at -fC < 1. 

For the initial state of the quantum map we take 
the Gaussian wave packet 



2N-1/4 



exp 



{x - XqY^ 
4(72 



iko{x - xq) 



where {x] 



4o^ 



((5a;2) EE - (.t)2 
and /cq is an integer. 



(9) 

= cr^, Po = {p) = 
The packet is as- 
sumed narrow: cr ^ and {5p^) ^ fikf). Note that in 
view of its periodicity in x the wave packet is generally 
not a state that minimizes the uncertainty relation. But 
in the case of a narrow packet it is essentially indistin- 
guishable from a minimum- uncertainty state |2^, ^ . 

A typical initial quantum state in studies of light 
squeezing is a coherent state [0, ^, |^. Such a state is 
an eigenfunction of the annihilation operator a, which in 
the present notation can be written as 



a = -^{y/^x + j-7=). 



2h 



1/2 



(10) 



we find it convenient to take the values of x and p immediately 
after the n-th kick rather than before the n-th kick as is done 
in majority of papers. It must be noted however that nonlinear 
dynamics is, of course, same for both these forms of the standard 
map. 
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The fact that the annihilation operator has such an 
appearance can easily be understood if we consider the 
following limiting case of the harmonic oscillator that 
follows from (H): 



if) = . 1 III 

dt 2 dx^^ 2 



(11) 



Now we can show that the wave function (Q) is a coherent 
state, i.e., an eigenfunction of operator (10), if we put 

9 h 

- 



27 



(12) 



Let us now turn to the problem of squeezing. In light 
squeezing experiments [0 , the observable quantity is the 
variance of the generalized quadrature operator 



(13) 



where d is the phase of the reference beam in the ho- 
modyne detecting scheme. In the particular cases where 
6* = and = 7r/2, Eq. (13) yields the following expres- 
sions for the generalized position and momentum opera- 
tors 



= a 



X^ = -i{a-a^), [Xi,X2] = 2i (14) 



with the uncertainty relation {5X1) {5X2) > 1, where 
averaging is done over an arbitrary quantum state and 
equality is achieved for a coherent state. The standard 
definition of quadrature squeezing is the condition ||l|, |^ 

min{{5Xl),{5Xl)) <1, (15) 

i.e., the variance of one of the quadrature components is 
smaller than for the coherent state. 

In a more general case we consider the variance (SXg) 
of the operator (O), and the state is assumed squeezed 
if the value of {SJQ) in this state for some value of 6 is 
smaller than in the coherent state . Experiments 

actually determine the minimum S of this variance as a 
function of the angle 9 as 



5= min (SXh 

9e[0,2Tr] 



(16) 



Using the definition ( |l^ ) of Xg, it is possible to show 
M M that 



and the minimum of (SXg) is reached at an optimum 
phase value 9 — 9* defined as follows [|l) : 



(18) 



For our discussion it is convenient to express S in terms 
of the cumulants of the operators x and p. Using the 
definition ([l0| ) of operator a and Eq. (|l7|), we obtain 



+ (fe')7 - v/(('5^'>7- (V)/7)'+4c2^ , 

(19) 



where c = ^{{{xp+px)) — 2{x){p)) Clearly for any Gaus- 
sian wave packet S — , while for a coherent state we 
have S = 1 because of equality (|l^). Hence a state is 
squeezed if 



S < 1 



(20) 



The condition determines the principal squeezing attain- 
able in homodyne detection j2^. The maximum of the 
variance (SXg) in 9 can be defined in the same way the 
minimum was defined in Eq. (p^. We denote it by S. 
Then we can show that the dependence of S on the cu- 
mulants differs from Eq. ( p^ ) only in the sign in front of 
the square root, so that we have 



SS> 1. 



(21) 



Thus, squeezing in S (Eq. po| ) is accompanied by dilation 
in S. Note that in contrast to the quadrature squeezing 
(|l5|), the definition of principal squeezing (19) contains 
quadrature correlators of the (xp) type. This is very 
important for systems with discrete time, to which the 
model of a quantum kicked rotator belongs. The thing 
is that the quadrature squeezing (15) is essentially unob- 
servable in such systems, although the principal squeez- 
ing ( p^ and ( pO| ) may occur ^. We will discuss the time 
dependence of S in Sec. 



III. THE NUMERICAL METHOD 

Several features of the numerical method must be men- 
tioned. The interval in x from to 2t: is partitioned 
into N segments Aa: = 2tt/N, and the wave function 
^'(x) is represented by a discrete sequence of values (col- 
umn vector I vp)) of length N, so that = ^'(^Aa;), 
I G [0,1, . . . , N —1]. Accordingly, k varies from to A^ — 1 
in the sum (^). In our numerical method N is always an 
integral power of two and the operator F in Eq. (j^) is 
represented in the form of the fast Fourier transform in- 
ducing the transformations 



(22) 



S^l + 2{6a^Sa)^2^{Sa^){Sa^^) (17) f^J 



To determine the principal squeezing, we must calculate 
((5(j2), {6p'^), and (xp) (see Eq. (19)). For instance, the 
calculation of (xp) proceeds along the following lines 

(xp) = (^\xF-^pF\^), 



where {'i>\ is obtained by transposing the vector \'i>) 
and then finding the complex conjugate of the result. 



^ Note that Lan |23| have studied the time dependence of {Sq"^) 
and (<5p^ ) for a quantum kicked rotator on a time scale on which 
one-to-ori£_nuantum-classical correspondence holds (see Table I 
in Ref. |23|). Both these variances increase uniformly so that 
quadrature squeezing llSh is impossible 
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while X and p are vectors that initially have the form 
x= [0,Ax,2Ax,...,2n- Ax],p^ [0, 1, 2, . . . , - 1]. 

The fact that x is defined in (mod 2tt) requires 
following the wave packet and ensuring that it is de- 
fined correctly during the passage through the end- 
points of the interval [0,27r]. We set up the process in 
the following manner. When the center of the wave 
packet in the ^-representation approaches an edge of 
the half-interval [0, 27r[, the wave function ^'(a;) is ex- 
amined on a new interval, [— 7r,7r[, with a new vector 
X = [0, Ax, . . . , TT, — 7r-|- Ax, — 7r-|-2Ax, . . . , —2 Ax, — Ax], 
since (— fcAx) mod 27r = (27r — kAx) mod 2tt, where k 
is an integer. The transition from [— tt, 7r[ to [0, 27r[ is 
treated similarly. 

Calculations in the p-reprcsentation have their own 
special features. Although for the Hamiltonian (|l|) the 
momentum is defined in the interval from — cx) to +oo, 
in numerical calculations we deal only with a finite range 
of values of momentum p, which is specified by the num- 
ber terms N in the Fourier transform (^. To avoid the 
possible problem of reflection of the wave packet from 
an edge of the given interval in the p-representation'^, we 
select this interval in each iteration of map (^) in such a 
way that the maximum of the absolute value of the wave 
function of the packet is always at the center of the given 
interval (actually, we renumber the vector p) . 

The process of calculating the next iteration of the 
quantum map (^) is terminated as soon as the wave 
packet ceases to be sufficiently localized either in the x- 
representation or in the p- representation, i.e., when the 
number of terms in the Fourier transformation actually 
involved in the calculation process is smaller than needed. 
Let us now to present the conditions for wave packet 
delocalization, which we use in our numerical simula- 
tions. First, introduce the notations ^ = max[o.27r] 
X = max{|Ai|, IA2I, . . . , |Ajv|}. Next, introduce Aiejt 
i-^right), thc valucs of Ak belonging, respectively, to the 
left and right edges of thc finite interval in which the 
wave function in momentum space is determined. The 
wave packet is delocalized and further calculations are 
terminated when one of the two inequalities, 

max{\Aieft\/x,\Arzght\/x} > £ or |*(z)|/^ > e 

become valid (here z = 0, if x € [0, 27r[ or z = tt if 
X G [— TT, 7r[). In this paper we used the cut off value 
£ = 0.002. 



IV. THE MAIN RESULTS 

As an initial wave function in our calculations we took 
the coherent state (a Gaussian wave packet) with h — 



^ For a discussion of the problem of reflection and splitting of a 
wave packet due to the finite range of momentum in the close 



10~^, ko = 10000, and a was varied between 0.04 and 
0.07. We fixed the initial width of the wave packet a and 
the value of Chirikov parameter K. Then, the parameters 
a and (3 in the evolution operator (||) are 



2cr2 



2a' 



(23) 



These formulas are obtained by combining the definition 
K :^af3 and Eqs. (0) and (|l|). In Sec. II we found that 
thc number of unperturbed levels of rotator captured in 
one kick is « 2a /h. From ( ^ ) it follows that in our case 
this number is K^l'^ j a' and amounts to several tens of 
thousands for the adopted widths of the wave packet. 

In our calculations K was varied between 0.2 and 2 
with a step of 0.02. We found the time dependence of the 
squeezing S ( p9| ) and the optimum value of the phase Q* 
at which {SXj) is at its minimum. To demonstrate the 
correlation that exists between the degree of squeezing 
and characteristics of chaos [|[ Q we calculated 



[{5x') + (V)] 



2\ll/2 



(24) 



It can be shown ^ that in the classical limit and 

while the wave packet is well-localized, i.e., [((5x^)]^/^ ^ 
xq and [(f^p^)]^^^ <C poi the d of ( p4[ ) corresponds to the 
following distance in phase space 



4/(t) = [{Axf + {Apn'\ 



(25) 



where (Ax(t), Ap(i)) is the solution of the linear small- 
perturbation equations near the classical trajectory 
{x{t),p{t)). The quantity dci{t) characterizes the diver- 
gence of two initially close trajectories and enters into 
the definition of the largest classical Lyapunov exponent 



A = lim 



dciit) 



(26) 



For a classical standard map in conditions of strong 
chaos, K ^ 1, there exists the simple dependence A « 
ln(X/2) Lyapunov exponent ( |2^ ) is an asymptotic 
characteristic of chaos. For finite time intervals 



dci{t) ~ exp(/i(x,p)t). 



(27) 



model of the quantum Arnold cat see Ref. p^ . 



where the exponent ft, is a function of a point in phase 
space and coincides, in order of magnitude, with the Lya- 
punov exponent A, but in some time intervals the differ- 
ence between the two may be significant. The latter fact 
can be explained by the strong inhomogeneity in the sta- 
tistical properties of the phase space of chaotic systems 
and, correspondingly, by the different rates of divergence 
of trajectories in different regions of phase space through 
which the system passes in its time evolution. It must 
be noted at this point that the dependence of h on the 
parameter K is rather complicated. What is important, 
however, is only the property of the strong (exponential) 
increase of dd ( |27| ) in the presence of chaos, a property 
often called local instability [Q . When the motion is regu- 
lar, the time dependence of dd is much weaker - it follows 
a power function Q . 
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FIG. 2: Upper part of the figure: Time dependence of the 
logarithm of squeezing In S. Lower half of the figure: Time 
dependence of Ind (see Eq. (^)). Everywhere a;o = tt and 
(7 = 0.006. 



On the other hand, it is h that determines the rate of 
phase volume deformation: the stronger the local insta- 
bility, the greater the deformation of phase volume in a 
given time interval. Since in our case quantum-classical 
correspondence and the concept of chaos are well-defined 
only in a very short time interval, while the wave packet 
remains localized, it is meaningful to consider the cor- 
relations existing between the time dependence of the 
squeezing and that of the quantity d (see (p^), which in 
the classical limit becomes dd (see (Eq)). 

Figure ^ depicts the time dependence of the logarithm 
of squeezing In S and Ind for different values of when 
the center of the wave packet is initially located at the 
point xq = tt, po = = 0.01. This initial condition 
is close to a hyperbolic point through which the chaotic 
layer passes even when K is small (see Fig. Figure |^ 
shows that the larger the squeezing (the smaller the value 
of S) the larger the local instability (the larger the values 
of Ind). This valid up to n « 4, when the packet spread 
becomes so large that purely quantum effects become im- 
portant. 

For another initial condition, xq = 7r/2, po = 0.01, 
which is closer to an elliptic point and hence lands in the 
chaotic region only at large values of if, the dynamics of 
squeezing is depicted in Fig. |[ We see that in this case 
squeezing is in about two orders stronger than in the 
conditions of Fig. ^ in the same time interval. On the 
other hand, both Fig. || and Fig. ^ exhibit an increase 
of squeezing as a function of the parameter K, which 
controls the strength of chaos in the system. 

Now we turn to study of the correlation between 



* We should note that for K = 0.8 (Fig. |l|) this position of wave 
packet does not belong to chaotic layer. 
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FIG. 3: Same as in Fig. ^ but for xq = tv/2. 

squeezing and the degree of local instability at differ- 
ent K in greater detail. The if-dependence of the de- 
gree of squeezing calculated after a fixed number of kicks 
at = TT and po = 0.01 is depicted in Fig. ^. After 
the third kick the correlation between In S and In d be- 
come very evident (Fig. ^a). However, small discrepan- 
cies in this dependence may appear as the number of 
kicks grows. Such discrepancies become evident, for in- 
stance, after the fourth kick for 1.1 < A' < 1.4 (Fig. |b). 
After five kicks, n = 5, the correlation between InS" and 
Ind is restored (Fig. ^). Note that this behavior pat- 
tern is quite typical. Hence, to establish the correlation 
between local instability and squeezing more clearly, a 
certain procedure of coarsening (averaging) these quan- 
tities in the given time interval is needed. In our study we 
determine the minimal squeezing, Smm, and the maximal 
dmax in a time interval during which the packet remains 
well localized for most values of K considered here. We 
found that there is a distinct correlation between Smm 
and dmax- the larger the value of dmax the smaller the 
value of Smin, and vice versa. An example of such a de- 
pendence is depicted in Fig. ||, where Smin and dmax were 
calculated after six kicks. Note that the diagrams do not 
go farther than K > 1.7 because after six kicks the wave 
packet becomes delocalized for K > 1.7 and calculating 
averages and local instability becomes meaningless. 

We also studied the dependence of the dynamics of 
squeezing on the initial width of the wave packet a. The 
results are depicted in Fig. ^. Clearly, the narrower the 
packet the stronger the squeezing achieved in a fixed time 
interval. This dependence arises because a narrow wave 
packet travels farther along its classical trajectory than a 
wide packet, so that it undergoes stronger deformations 
related to nonlinear classical dynamics. The exponential 
decrease is replaced by growth when the wave packet de- 
parts from the classical limit and the dynamics is of an 
essentially quantum nature. 

Now examine the problem of stability and observabil- 
ity of squeezing in conditions of chaos. The figures men- 
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FIG. 4: The dependencies of InS (solid curves) and Ind 
(dashed curves) on the Chirikov parameter if at a fixed num- 
ber of kiclis: (a) n — 3, (b) n = 4, (c) n — 5. Everywhere 
a;o = TT and a = 0.007. 



tioned earlier can serve to illustrate the statement that 
the stronger the chaos the stronger the principal squeez- 
ing. However, the definition of principal squeezing (W|) is 
related to fixing the phase, 9 ~ 6* . Here 9* is time depen- 
dent even for exactly integrable systems ||2l| . For strong 
chaos in the classical limit, the time dependence of 9* 
may be very complicated. Indeed, in addition to stretch- 
ing and squeezing, the main feature of classical chaos in 
the systems with bounded phase space is the multiple for- 
mation of folds of the phase volume as a system evolves 
0. Hence the process of finding the "minimum width" 



FIG. 6: Time dependence of the logarithm of squeezing S 
for different initial widths of the wave packet a and at fixed 
values K — 0.8, xo = tt. 



of a phase drop, which actually amounts to finding the 
dependence 6* vs t in the quasiclassical limit, becomes 
unstable for large time intervals. 

Basing our reasoning on a similar semiclassical picture, 
we examined the stability of the time dependence of the 
optimal phase 9*{t), which was calculated quantum me- 
chanically, against a small perturbation of the initial po- 
sition of the wave packet. More precisely, we found the 
time dependence of the optimal phase 0* with the initial 
condition xq ^ n and, similarly, ^2 (^) with the initial 
condition xq = ir — 0.05. We denote the difference of 
these phases by 

Dit) = 9lit)-9;it). (28) 
Since 9*{t) is periodic with a period tt (see Eq. (Esl), it 
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FIG. 7: The difFerence of optimal phases D (see Eq. (|2^)) as 
a function of the Chirikov parameter A' at a fixed number 
of kicks: (a) n = 2, (b) n = 3, and (c) n — 4. Everywhere 
a;o = TT and a = 0.006. 

is natural to take sin(2Z?) as the quantity of interest, 
since m this way we avoid breaks in the diagrams related 
to the periodicity of 9*{t). The dependence of sin(2I?) 
on the Chirikov parameter K for different fixed num- 
bers of kicks is depicted in Figs. 0a-^. After two kicks 
(Fig. I^a) the maximum value of | sin(2_D) | does not ex- 
ceed 0.035 at K = 2. After three kicks ((Fig. |a) the 
value of D becomes significant at i^T > 1.2. Finally, af- 
ter four kicks (Fig. ^) the process of measuring squeez- 
ing becomes essentially unstable at X > 1. Indeed, in 
these conditions with a small perturbation of the initial 
position of the wave packet, the difference of the opti- 
mum phases reaches a value of order tt. Such regime of 
squeezed states generation was called unstable squeezing 
in As Fig. |^ implies, unstable squeezing is observed 
when chaos is strong. On the other hand, for short time 
intervals and small values of K, the squeezing can be 
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FIG. 8: Dependence of absolute value of the wave function 
I \1/ I on a; and dependence of the absolute value of the Fourier 
transform | Ak \ on k (see Eq. (^)) at (a) n = 6, (b) n = 18. 
Everywhere K = 1.2, xq = tt and a — 0.006. 

already strong enough but yet rather stable. 

To conclude this section we will briefly touch on the 
problem of the dynamics of disintegration of coherent 
states in conditions of chaos, a problem that is of inter- 
est by itself. Figures |8^ and ||b depict the dependence 
of I I on a; and of \ A \ on k (see Eq. (||)). Actually, 
Fig. H gives the shape of the wave function in the co- 
ordinate and momentum representations for an initially 
narrow wave packet with [{Sx^)]^^'^{t — 0) = a — 0.006 
and [(V)]i/2(i = 0) = (1/12) X IQ-^. The relatively 
small value K ^ 1.2 makes it possible to examine the 
fairly long evolution of the wave packet up to the point 
of its total disintegration^. After six kicks (Fig. ||a) the 
wave packet spreads out significantly, but on the whole 



^ In conducting numerical experiments in the dynamics of the dis- 
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retains its bell shaped structure. What follows is a disin- 
tegration of the packet into many small subpackets, with 
the characteristic shape of the wave function depicted 
in Fig. (after 18 kicks). Finally, very soon the wave 
function becomes so dissected that even 2^^ Fourier har- 
monics are insufficient to describe the evolution correctly 
(for the data of Fig. ^ this happens approximately at the 
20th kick). 

Qualitatively, the same pattern of the evolution of the 
wave packet was observed at higher values of K: first the 
broadening of the wave packet, and then its rapid disinte- 
gration into many very small subpackets. The differences 
in packet disintegration for large values of K in compar- 
ison with the case K ~1 (as in Fig. ^ boil down to two 
facts: (i) the "swelling of the packet" and the disintegra- 
tion occur very rapidly (it takes only several kicks to com- 
plete the process), and (ii) the emerging subpackets are 
extremely small. Hence the process of disintegration of 
wave packets in strong chaos resembles an explosion. On 
the whole, the pattern being described agrees well with 
the pattern obtained from the analysis of the behavior 
of the Wigner function |l8| , although we observed some 
anomalies. In particular, for fairly narrow wave packets 
(cr = 4 X 10~^) we observed the disintegration of the ini- 
tial packet into two fairly large subpackets. Ripples then 
appeared on the subpackets, and the two disintegrated 
into many small packets. A more detailed description of 
the disintegration of coherent states at chaos will be a 
subject of our separate publication. 



V. DISCUSSION AND CONCLUSION 

Thus, in this work we numerically study the dynam- 
ics of generation of squeezed states in the evolution of a 
Gaussian wave packet in the quasiclassical limit for the 
model of a kicked quantum rotator. We show that within 
the time interval where the packet is well-localized the 
squeezing becomes stronger in the transition to chaos. 
For strong chaos and in long time intervals, the squeez- 
ing process becomes unstable. These results, obtained 
through direct numerical simulation, are in good agree- 
ment with the results derived using perturbation theory 
for other models § § 0. 

In the final stages of preparing the manuscript for press 
we became acquainted with two recent papers [E9| also 

integration of wave packets we did not use the proceoure (de- 
scribed in Sec. ^ ) of terminating the counting process when the 



devoted to the problem of generating nonclassical states 
(squeezing and antibunching) in the systems with quan- 
tum chaos. The authors of Ref. |2^ presented the results 
of numerical experiments on the dynamics of quadra- 
ture squeezing in simple quantum models that allow a 
transition to chaos in the classical limit: the Lipkin- 
Meshkov-Glick model ]3C|] and the Belobrov- Zaslavsky- 
Tartakovskil model piT^ In contrast to our approach, 
the authors of paper |23| were interested in the long-time 
limit, when the wave packets are delocalized and this 
sense the quantum-classical correspondence is completely 
violated. They found that quadrature squeezing disap- 
pears at the transition to quantum chaos, although to 
some degree squeezing is always present in conditions of 
regular motion. The nonzero squeezing in the conditions 
of quantum chaos and in short-time limit has been also 
mentioned in |^9| , but the enhance of squeezing was not 
observed, probably because in their numerical simula- 
tions the quasiclassical parameter was not large enough: 
only several hundred of quantum levels participated in 
the dynamics of the system. Thus, the results of |2|] 
do not contradict ours but supplement them in another 
limiting case, the limit of long times of motion. 

In conclusion wc would like to make a remark concern- 
ing the possibility of observation of squeezing in quan- 
tum chaos on a time scale corresponding to a well-defined 
quantum-classical correspondence. At present essentially 
all experiments on light squeezing are done in the sta- 
tionary regime. Squeezing at the transition to quantum 
chaos increases only over finite time intervals and in this 
sense is a transient dynamical phenomenon. We hope 
that the development of effective experimental methods 
for observing squeezed states of light in transient dynam- 
ical regimes will also make it possible to observe the en- 
hanced squeezing at the transition to quantum chaos. On 
the other hand, as noted in the Introduction, it is much 
simpler to realized the quantum kicked rotator model in 
atomic optics |27|. Moreover, it is much simpler to ob- 
serve transient dynamical regimes in experiments with 
cold atoms. Hence we believe that atomic optical sys- 
tems have great potential for observations of squeezed 
states at transition to quantum chaos. 

We are grateful to Audrey Kolovsky and Jan Pefina for 
discussions and to Boris Chinkov for support and drawing 



our attention to his work 1 10 



wave function becomes delocalized. 
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